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O Abstract 

We present results on modeling solar pores and sunspots using 2D axisymmetric 
' magneto-hydrodynamic (MHD) simulations. These models are helpful for understand- 
ing the mechanisms of magnetic field concentration in sunspots, and the large-scale 
flow patterns associated with them. The simulations provide consistent, self-maintained, 
^ although not fully realistic, models of concentrated magnetic field near the solar surface. 
\^ In this paper, we explore under which conditions the associated flows are converging 
l/^ I or diverging near the surface. We find that in most cases in which a stable, pore-like 
concentration of magnetic field forms, a configuration with converging over diverging 
flow is established. 
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1. Motivation and Objectives 

The mechanisms of how magnetic pores and sunspots form on the Sun are still poorly 
understood. Observations and numerical simulations suggest that their structure is inti- 
mately linked with characteristic surface and subsurface flows in and around the mag- 
netic field concentration. Fully developed sunspots exhibit a surface outflow in their 
penumbra, the so-called Evershed flow, presumably caused by interaction between the 
near-surface granular convection and the highly inclined penumbral magnetic field as 
suggested by the numerical simulations of Kitiashvili et al. (2009). Observation (e.g., 
Zhao, Kosovichev, and Sekii, 2010) have revealed downflows in the central region of 
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the sunspot and subsurface converging flows (inflows) below the granulation layer, as 
well as outflows further below. Numerical models of subsurface magneto-convection 
in axisymmetric configuration (Hurlburt and Rucklidge, 2000; Botha, Rucklidge, and 
Hurlburt, 2006; Botha et ai, 2008) have been able to reproduce similar flow struc- 
tures. In most of their simulations, inflows over outflows were found. It is conjectured 
that this flow is fundamentally important for maintaining the integrity of the mag- 
netic field concentration. Recently, a numerical simulation of realistic compressible 
magneto-hydrodynamic (MHD) convection has been reported in which a magnetic field 
concentration was held in place by boundary conditions (Rempel et ai, 2009). A real- 
istic Evershed flow was found at the surface, but the outflow continued well below the 
surface. Botha, Rucklidge, and Hurlburt, 2006 and Botha et ai, 2008 did in fact also 
found cases with a diverging flow over a converging flow. The present paper extends 
their work exploring in more detail under which conditions diverging or converging 
flows can hold a magnetic flux concentration in place. 

2. Numerical method 

We study magneto-convection in an axisymmetric cylindrical geometry using a code 
originally developed by Hurlburt and Rucklidge (2000) for the two-dimensional (2D) 
case and later extended by Botha et al. (2008). The exact equations and their parameters 
are discussed in detail in these two publications. The model considers a layer of electri- 
cally conducting, perfect monoatomic gas subject to uniform gravitational acceleration, 
with constant shear viscosity, magnetic diffusivity, and magnetic permeability. It ap- 
proximates the conditions in the upper part of the solar convection zone but excluding 
the very top few hundred kilometers below the photosphere where the plasma is only 
partially ionized and radiation needs to be accurately modeled. For thermal boundary 
conditions, we prescribe a constant heat flux at the bottom and Stefans law at the top. 
The side wall is perfectly electrically conducting and does not allow for a heat flux 
across it. Top, bottom, and outside walls are impenetrable and stress free. The magnetic 
field is vertical at the bottom and matched to a potential field at the top. 

The computational domain can be subjected to constant rotation by including Cori- 
olis and centrifugal forces. The equations are written and solved in non-dimensional 
form with the key parameters being: Rayleigh number, R, Prandtl number, a, magnetic 
Prandtl number, CJ„„ temperature contrast between the top and bottom of the domain, 
0, rotation rate, H, aspect ratio, the ratio between height and radius of the cylindrical 
domain, F, and the Chandrasekhar number, Q, a measure of the magnetic flux in the 
system. The ratio between specific heats is chosen to be y = 5/3, appropriate for a 
monoatomic ideal gas. 

The initial temperature and density profiles in the simulations take the form of a 
polytrope, in non-dimensional form T{z) = I + 6z,p{z) = (1 + 9z)"\ where T, p, z, 
and m are the non-dimensional temperature, density, depth (ranging from at the top 
of the domain to 1 at the bottom), and the polytropic index, respectively. Simulations 
are started with an initial uniform vertical magnetic field and are run until the system 
reaches a more or less steady state, if such a state is attained at all for the given set of 
parameters. The results presented in this paper are all for an aspect ratio of F = 3. 

The equations are solved using a finite-difference scheme accurate to sixth order and 
a fourth-order time marching scheme. 
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Figure 1. Flow visualization (arrows) and strength of the vertical magnetic field (contours) for a case that 
forms a diverging flow over a converging flow. See text for parameters. 
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Figure 2. Dependence on the rotation rate. Results are shown for the same parameters as the reference 
case (figure 1 ) except the rotation rate, £2, is reduced here from 0. 1 to 0.02 (left panel) and (light panel), 
respectively. 



3. Results 

This work is an extension of Botha, RuckHdge, and Hurlburt, 2006 and Botha et ai, 
2008 in which we are exploring the conditions under which a stable flux concentration 
with an outflow over an inflow can exist. Unless otherwise specified, the simulation 
parameters for the results presented here are: Q = 32, m = 1, F = 3, = 0.2, i2 = 0.1, 
(7 = 1, = 10, and R = 10^, referred to as the reference case in the text below. This 
is one of the few cases in which a diverging over converging flow structure forms. 
Simulation results for the same parameters were originally presented in figure 17 of 
Botha, Rucklidge, and Hurlburt, 2006. A visualization of the flow and the magnetic 
field strength in shown here in figure 1. It seems that externally applied rotation is 
crucially important in maintaining this flow structure. Figure 2 shows two examples in 
which we reduced the rotation rate resulting in a weakening of this flows configuration. 
In fact, we have not found a parameter set without rotation in which a stable outflow 
over inflows was established that at the same time tightly held the flux concentration 
together. 

Similarly, reducing the Prandtl number, a, to more realistic values weakens this 
outflow configuration and a configuration with inflows over outflows, more similar to 
the observational evidence, is found. This is illustrated in figure 3. 
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Figure 3. Prandtl number dependence. Results are shown for the same parameters as the reference case 
(figure 1) except the Prandtl number, a, has been reduced from 1.0 to 0.3 (left panel) and 0.1 (right panel), 
respectively. 
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Figure 4. Same parameters as the reference case except for 9 = 20 and i? = 4 x 10^. 



Lastly, increasing the strength of the convection by increasing the Rayleigh number, 
R, also causes the flow to establish the more familiar converging over diverging flow as 
shown in figure 4. A weak secondary convective cell is formed in this case further away 
from the axis. 



4. Conclusions 

Admittedly, the simulations presented here are somewhat limited because of their sim- 
plifications and the restricted geometry, but they nevertheless can give us some insight 
into the mechanisms that are involved in maintaining a tightly concentrated magnetic 
field near the solar surface inside structures such as magnetic pores or sunspots. In 
fact, due to the simplicity of these simulations, a whole range of parameters can be 
explored rather easily. The simulations have shown that for cases in which a stable, 
tightly concentrated magnetic structure is formed, most often a flow configuration with 
an inflow (towards the spot axis) above an outflow deeper below is established, compat- 
ible with observational evidence. It is conceivable that such a configuration is better for 
maintaining the magnetic flux concentration at the center since the dynamic pressure of 
the inflow helps balancing the magnetic pressure inside the flux concentration. A flow 
with opposite sign is not impossible, though. In fact, under certain conditions a outflow 
configuration is estabUshed. Rotation seems to be crucially important in these cases. 
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